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EFFECT  OF  UNEgUILIBRIZED  DISSOCIATION  ON  FLOW 
AROUND  BLUNT-NOSED  BODIES 
by 

Yu.  P.  Lun'kin  and  M.  P.  Shtangel' 


In  this  study  we  consider  the  effect  of  unequilibrixed  dissocia¬ 
tion  of  a  diatomic  gas  on  the  position  of  the  shock  wave  and  the  pro¬ 
file  of  the  gas-dynamic  pjarameters  along  the  zero  line  of  the  flow  of 
a  blunt-nosed  body. 

1.  Basic  Equations 

Let  us  consider  the  flow  around  an  axisymmetrical  body  with  a  sphertcal  head 
part  by  a  diatomic  disassociating  gas.  Since  the  time  of  excitation  of  the  revol¬ 
ving  and  oscillating  degrees  of  freedom  is  much  less  than  the  time  of  the  dissocia- 
•  ion  relaxation,  therefore  in  agreement  with  the  method  of  quasi-equillbrized  zones 
Ll]  it  is  possible  to  consider  that  the  inner  degrees  of  freedom  in  the  process  of 
dissociation  find  themselves  in  a  state  of  equilibrium. 


Here  v  and  Vr  are  components  of  the  speed j 

P»  density,  pressure,  temperature  and  degree  of  dissociation; 

D  is  the  energy  of  dissociation; 

C,  are  constants  which  characterize  the  respective  speed  of 

dissociation  and  the  properties  of  the  concrete  gas. 

The  equation  of  the  energy  along  the  line  of  the  flow  and  the  equation  of  the 
state  will  have  the  form 
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where ^  is  the  molecular  weight  of  the  diatomic  gas; 

i^  and  i^  are  the  molar  enthalpies  respectively  of  the  monoatomic  and  diatomic  gases 
and  i^^  is  the  enthalpy  of  retardation. 

Let  us  note  th^t  when  there  occurs  dissociation  of  the  molecules,  the  oscilla¬ 
tions  are  already  considerably  excited  and  the  heat  capacity  of  the  molecules  prac¬ 
tically  does  not  depend  on  the  temperature,  i.  e.,  i^  =  Cp2T  and  i^^  =  CplT. 

Let  us  introduce  the  new  variables 
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Then  the  system  (1)  -  (5)  in  the  new  variables  will  be  written  (the  line  over  the 
values  will  be  omitted  by  us) 
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Here  by  the  exponent  2  there  are  designated  the  values  after  the  shock  wave. 

p  ^  (f ,,  p.) 
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where  'Jila  the  parameter  of  the  inequality  withij"^!  1  flow  pr''ctically  frozen 
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and  with  ^  i  flow  close  to  equilibrium. 

One  should  note  that  the  equations  (B)  -  (iO)  remain  unchangeable  with  the  sub¬ 
stitution  of  X  and  ax,  and  the  equation  (L2)  does  not.  This  points  to  the  fact  that 
in  the  equilitrium  case  the  value  for  the  withdrawal  of  the  shock  wave  is  propor¬ 
tional  to  Rq  but  the  unequilibrized  flow  does  not  have  this  similarity.  For  the  os¬ 
cillation  relaxation  this  was  noted  in  the  reference  [2]. 

2.  Method  of  Solution 

We  will  seek  a  solution  of  the  system  (8)  -  (12)  in  the  form  of  series  by  degree 
-tf'with  coefficients  depending  on  x  [2],  [3]* 

.7  II,.  (  \  1  H  U  (A  )  ... 

V  i',  (.V)  •  1'  (A  )  ^'  t'  •  •  • 

P  ^  (  A  )  r  P  ( V)  t  ... 

P  P„t-v)  •  p  (A  )  -r  •  • 

T  7„(a  i  r-  T  (  V)  I-  .  . . 

1  1  (a  )  •  1  ^a)  . 

By  substituting  the  expressions  (15)  in  the  equations  (8)  -  (12),  we  shall  get 
the  infinite  system  of  ordinary  differential  equations  for  the  detenniiiation  of  Uj^, 

Vfi  ,  ar.a  30 

In  the  lowest  approximation  in  accordance  with  -O' we  have  (the  subscript  0  we 
vvlll  orrlt) 
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Aith  a  sufficient  degree  of  precision  one  can  accept  that  close  to  the  zero  line  of 
flow  the  shock  wave  proves  tc  be  a  sphere  concentric  v'lth  the  body  [3].  L«t  us  in¬ 
troduce  the  dimensionless  radius  of  the  shock  wave  _ i  i  ^  ^ 
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With  X  -  y\  we  have  in  accordance  vdth  [i],  [3] 
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where 
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The  systeiTi  (id)  -  (20),  (13)  does  not  turn  out  to  be  closed  since  it  contains 
seven  unknowns.  Numerical  computations  showed  that  both  for  the  equllibrized  and 
for  the  nonequilibrized  flow  v(x)  on  the  zero  line  of  the  flow  changes  almost  lin¬ 
early  iu).  Therefore  by  stating  v(x)  in  the  form, 

■  u' (/.'I  (23) 

we  will  have  a  closed  system  of  equations  for  the  determination  of  v,  u,  p,  pi,C(  , 

Pi¬ 
ny  malcing  use  of  the  system^  (16)  -  (20),  (13)  and  the  boundary  conditions  (21) 
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If  in  the  equation  (23)  we  substitute  the  relationships  (24),  (25)  and  put  x  =  1 
(v(l)  =  0),  then  we  will  get  an  algebraic  equation  of  the  fourth  degree  relative  to 
Hence  it  is  easy  to  fina  the  value  for  the  withdrawal  of  the  shock 
wave  along  the  zero  line  of  flow. 
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In  continuation  let  us  make  use  of  the  fsct  that  the  term 
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2  in  the  equation 

(^3)  is  small  in  comparison  v<rith  the  retrAining  terms  [5],  and  it  can  be  disregarded. 

Besides  this  the  pressure  along  the  zero  line  changes  insignificantly  and  it  can  be 

3  P 

considered  as  constant  and  equal  to  p  =  1  -  ^  in  accordance  with  [6]  Then 

from  the  equations  (IV),  (20),  (13)  we  shall  find  =  C)f(x),  T  =  T  (x)  and /^  =  /^(x). 


In  particular  for  we  have  the  expression 
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Afterwards  from  the  equation  (17)  it  is  possible  to  determine  the  profile  of 

the  pi-essure  and  find  the  second  approximation  fcr  flf  (x),  f  (x),  T(x) 

c'l . . --ady  with  the  taking  into  account  of  the  corrected  pressure.  However,  the  expres- 

/ 


sions  obtained  prove  to  be  quite  cumbrous. 


I'igure  1. 


Figure  2. 
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Figure  3. 
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3.  Evaluation  of  the  Results  Obtained 

A  program  was  compiled  for  the  BESM-2  and  computations  were  made  for  a  wide 
range  of  parameters  of  the  Inequality  of  ^and  different  figures  for  M,  From  the  re¬ 
sults  obtained  one  can  draw  the  following  basic  conclusions: 

1.  The  value  for  the  departure  monotonlcally  decreases  with  the  Increase  In 
the  M  number  with  a  given  CRq  (C  is  the  constant  whilch  characterizes  the  speed  of  the 
reaction,  and  R^  is  the  radius  of  the  curvature  at  the  critical  point)  and  Increases 
with  the  decrease  in  CR^  with  a  given  M. 

2.  For  the  interval  under  consideration  of  the  M  numbers  and  the  parameter  of 

Inequality  profiles  of  speed  prove  to  be  practically  linear. 

3.  The  profiles  of  temperature,  density  and  degree  of  dio:;ociation  substantial¬ 

ly  depend  on  X-  wit-h  Ixj  ^  1  the  basic  change  in  these  values  occurs  close  to 
the  body,  and  with  1  right  behind  the  shock  wave. 

As  an  examjjle  In  Figures  1,  2  and  3>  respectively,  profiles  of  density,  temper¬ 
ature  and  degree  of  dissociation  for  M  number  20  at  CR^  =  U  *  lo”^,  !~X( ~  2.1683  •  10 
CRq  -  4  •  10^,  IXh  0*21683,  CRq  «  4  •  10^°,  |X/=  2.1683  are  shown. 

In  conclusion  we  consiaer  it  our  duty  to  express  our  thanks  to  T.  Ya.  Timofeyev 
for  his  aid  in  compiling  the  program  and  carrying  out  the  computations. 
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